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We consider the dynamics of fermions with a spatially varying mass which 
couple to bosons through a Yukawa interaction term and perform a consistent 
weak coupling truncation of the relevant kinetic equations. We then use a 
gradient expansion and derive the CP-violating source in the collision term 
for fermions which appears at first order in gradients. The collisional sources 
together with the semiclassical force constitute the CP-violating sources rele- 
vant for baryogenesis at the electroweak scale. We discuss also the absence of 
sources at first order in gradients in the scalar equation, and the limitations 
of the relaxation time approximation. 
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1 Introduction 



The main unsolved problem of electroweak baryogenesis JI| is a systematic computation of 
the relevant sources in transport equations. We shall now present a method for controlled 
derivation of leading CP-violating sources appearing as a consequence of collisions of 
chiral fermions with scalar particles in presence of a scalar field condensate. We assume 
the following picture Q] of baryogenesis at a first order electroweak phase transition: 
when the Universe supercools, the bubbles of the Higgs phase nucleate and grow into 
the sea of the hot phase. For species that couple to the Higgs condensate in a CP- 
violating manner that CP-violating currents are created at the phase boundary (bubble 
wall). These currents then bias baryon number violating interactions mostly in the hot 
(symmetric) phase, where the B-violating processes are unsuppressed. The baryons then 
diffuse to the Higgs phase, where the B-violating interactions are suppressed, resulting in 
baryogenesis. 

Kimmo Kainulainen || has explained how to systematically derive the CP-violating 
source in the flow term of the kinetic equation for fermions. For details see Paper I ||. 
The source is universal in that its form is independent on interactions. It can be rep- 
resented as the semiclassical force originally introduced for baryogenesis in two-Higgs 
doublet models in 01], and subsequently adapted to the Minimal Supersymmetric Stan- 
dard Model (MSSM) in ||. This problem involves computation of CP- violating sources 
from charginos, which couple to the Higgs condensate in a manner that involves fermionic 
mixing. Here we show how to compute the CP-violating source in the collision term 
that arises at first order in gradients. We work in the simple model of chiral fermions 
coupled to a complex scalar field via the Yukawa interaction with the Lagrangian of the 
form H 11 1 

C = i^ppip - ip L mip R - ip R m*ip L + £ yu , (1) 
where £ yu denotes the Yukawa interaction term 

£ yu = -vW>i$r - y<P*^R4>L, (2) 

and m is a complex, spatially varying mass term 

m{u) = y'§ = m R (u) + im^u) = \m(u)\e l9( - u \ (3) 

Such a mass term arises naturally from an interaction with a scalar field condensate 
$o = This situation is realised for example by the Higgs field condensate of a first 
order electroweak phase transition in supersymmetric models. When (f> in @) is the Higgs 
field the coupling constants y and y' are identical; our considerations are however not 
limited to this case. 

The dynamics of quantum fields can be studied by considering the equations of motion 
arising from the two-particle irreducible effective action (2PI) |7| in the Schwinger-Keldysh 
closed-time-path formalism || [5| . This formalism is for example appropriate for studying 
thermalization in quantum field theory 
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Figure 1: The Schwinger closed-time-path (CTP) used in the derivation of the 2PI 
effective action (g). 

We are interested in the dynamics of the fermionic and bosonic two-point functions 

S a p{u,v) = -i (T c [ipa{u)ip(3{v)}) (4) 

A(u,v) = -i(T c {cf>{u)<t>\v)]) (5) 

where the time ordering Tq is along the Schwinger contour shown in figure [I], which is 
suitable for the dynamics of out-of-equilibrium quantum fields. 

2 Effective action and self-energies 

The 2PI effective action can be in general written in the form |7]. || 

r eff = r + Fx + r 2P i, (6) 

where the tree-level and one-loop actions read 

T = i [ d%d 4 v[A 1 (u,v)A(v,u) - S 1 (u,v)S(v,u)} (7) 
Jc 

r\ = t J c d%[\n A- 1 {u,u) -In S-^u^)}, (8) 

and Aq 1 and Sq 1 are the tree-level inverse propagators: 

S \u,v) = {i$ u -m{u)P R -m*{u)P L )5 c {u-v) (9) 
Ao 1 ^) = (-dl-ml(u))5 e (u-v). (10) 

Here Pl,r = (l = F7 5 )/2, and 5c is the ^-function along the contour C. F 2 pi in (^) contains 
higher order quantum corrections, which can be for example studied in the loop expansion. 
At two loops F 2 pi reads 

r 2P i - F 2 [S, A] = -y 2 J c d 4 ud 4 v Tr[P R S(u,v)P L S{v,u)]A(u,v), (11) 

which can be easily computed from the two-loop diagram in figure ||. The one-loop self- 
energies are then obtained from ( pTTj ) by taking the functional derivatives 

V(u,v) = -i /J™ ; = iy 2 [PlS{u, v)P r A{v, u) + P R S{u, v)P L A(u, v)] (12) 
n'(u,v) = * 5 ^ l u) = -ty 2 Tr[P L S(v,u)P R S(u,v)}. (13) 




Figure 2: The diagrams contributing to the 2PI effective action (|7HIT|) up to two loops for the 
Lagrangian (|IJ) with the fermion-scalar Yukawa coupling term (|2|). The full (dressed) fermionic 
and bosonic propagators are denoted by S (solid blue lines) and A (dashed red lines), respectively. 
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Figure 3: The Dyson- Schwinger equations at one loop obtained by varying the 2PI effective 
action (|7nl~l]) with respect to the fermionic and bosonic propagators. The corresponding tree- 
level Lagrangian is given by The full (dressed) fermionic and bosonic propagators are 
denoted by S (solid blue lines) and A (dashed red lines), respectively. 

The equation of motion for S and A are obtained by varying the effective action with 
respect to S and A. The resulting equations are simply 

S^^S = 5 C + Z'®S (14) 
A^fgiA = 5 C + IT® A, (15) 

where ® denotes a convolution with respect to contour integration. These are the Dyson- 
Schwinger integro-differential equations diagrammatically shown in figure []. The simple 
look of these equations is deceptive, since they involve integration over the closed-time- 
path in figure [I]. To proceed we use the Keldysh reformulation of the problem, according 
to which the contour C is split into two parts: (to,t) in the positive time direction and 
(t, to) in the negative time direction; finally we set to ~~ * —oo. This corresponds to the 
replacements 
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S{u,v) — > S ab (u,v) 
A(u,v) — ► A ab {u,v) 
S c (u-v) — > ad ab 8(u-v) (16) 

in the effective action (|BT-|i~T|). This procedure naturally leads to the Keldysh 2x2 formu- 
lation for the two-point functions (f|-|5|) in which the off-diagonal elements of the fermionic 
Wigner functions correspond to 

S < (u,v) = S + ~(u, v) = i(^(v )i/j(u)) 

S > (u,v) = S- + (u,v) = -i(ip(u)$(v)) } (17) 
and similarly the bosonic ones @ are 

A < (u,v) = A + -(u,v) = -i(<p ] (v)<p(u)) 

A>(u,v) = A- + (u,v) = -i(<f>(u)<fi(v)). (18) 

In the Keldysh representation it is convenient to redefine the self-energies ( |i2| - |l"3"|) as 

Z ac {u,v) = aZ' ac {u,v)c (19) 
U ac (u,v) = aW ac (u,v)c. (20) 

3 Kinetic equations 

When written explicitly the kinetic equations for fermions fll4]) and bosons ( pT5[ ) become 
^ + ^-mP R -m*P L )S < -e- i<> {J: R ,S < }-e- l<> {^ < ,S R } = (21) 
(-^ 2 + A; 2 + ^-,9-m 2 )A < -e-^{n i? ,A < }-e-^{n < ,A R } = C+, (22) 



S K {x,k) = I d 4 re lk - r S < (x + r/2,x-r/2) (23) 



where 



defines the Wigner functions for fermions, and there is a similar expression for bosons. IIr, 
Sr and A R denote the hermitean parts of the self-energies and two-point functions, 
respectively. We postpone a discussion of the physical relevance of the contributions 
involving the hermitean parts to a later publication. Eq. (|23| ) defines a Wigner transform, 
which is useful to separate the dependence on the slowly varying average (macroscopic) 
coordinate x — (u + v)/2 from that on the microscopic coordinate r = u — v. 
The bilinear operator o in Eqs. ( |2l"l - p2|) is defined as 



o {A, B}= 1 - (d x A ■ d k B - d k A ■ d x B) , (24) 
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the mass terms read 



rh = m(x)e i d ' 9k (25) 

ml = m\{x)e-^ d ' d \ (26) 

and the collision terms are of the form 

= -I C -*({E> S<} (27) 

Q = -V^({n>,A<}-{n<,A>}). (28) 



The one-loop expressions for the self-energies can be inferred from Eqs. ( |iT?H2"0| ) and (|T2|- 



13): 



£<'>(£, x) = E + -~ + (k,x) 

+ {2v)%k — k' — F)P/ ? 5 < ' > (A; / , s)P i A<' > (A;", x)] (29) 

!!<'>(*;, x) = n+"'" + (A;,x) 

= ~iy 2 J ^^(2n)H(k + k'-k'')Tr[P R S>'<(k',x)P L S<>(k'',x)}. (30) 

4 Wigner functions 

We are interested in modeling the dynamics of fermions and bosons in the presence of a 
Higgs field condensate of growing bubbles at a first order electroweak phase transition. 
When the bubble walls are thick we can expand in gradients of the condensate. This 
expansion is accurate for quasiparticles whose de Broglie wavelength E^b is small when 
compared to the scale of variation of the background, which is specified by the phase 
boundary thickness L w . Since equilibrium considerations yield L w ~ 5 — 15, and the de 
Broglie wavelength is typically given by the inverse temperature, Ids ~ we have 

^dB^x ~ £<ib/L w 1, so that the expansion in gradients is justified. 

The problem can be further simplified by noting that typically large bubbles are almost 
planar and slow, such that it suffices to keep the leading order terms in the bubble wall 
velocity, and expand to leading nontrivial order in gradients. So when written in the rest 
frame of the bubble wall (wall frame), Eqs. ( |2"5j - Po] ) simplify to m(z) = m + \m!dk z + .. 
and ml(z) =ml + ^m^d kz + .. . 



In Paper I we have shown that the fermionic Wigner function (T?) acquires a nontrivial 
contribution at first order in gradients. In the presence of a scalar condensate, when one 
neglects particle interactions, a wall moving in ^-direction conserves spin in z-direction. 
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This implies that one can write (cf. Paper II) the Wigner function in the following 
spin-diagonal form 

s<>> = J2 s; 



r<,> 



s=± 



where 



S<>> = iP s 



3 5 s<,> 3 s<,> . s<,> • 5 s<,> 

si 7 9o -7% +9i ~^92 



Here P s denotes the spin projector 



(31) 
(32) 



'l + sS z 



(33) 



where 



S-, 



r (A;o7 -£|| -7||)7 3 7 5 
fc 



(34) 



is the spin operator in z-direction, and k$ = sign[fco]y &o — kf\- 

The quantity ' in Eq. (|32|) is a measure for the phase space density of particles 
with spin s. When written in the wall frame and to first order in gradients, the solution 
for ^o <!> nas the form 



2ird(k 2 - \m\ 2 + ^s\m\ 2 6') \k \n{k 
v kn ' 



x 



St 



-2n6(k? 
In thermal equilibrium we have 

n(k, x) — > n = 



777, 



e k /T + 1 ' 



+ ^s\m\ 2 9')\k \(l - n(k,x)). 



k = 7 w {k + v w k z ), 



(35) 



(36) 



where v w = v w z is the wall velocity and 7^ = 1/yl — v^/c 2 . The quantities g 
(z = 1, 2, 3) in (HID can be related to Qq ^ by (see Paper I & II) 



s<,> 
i 



s<,> 

9i 



s<,> 
92 



9^ 



1 

h 
1 

k L 
1 



s<,>\ 



mug^ - ■^d z {m I gl <>> ) - S -^ L d kz (k z g s 
2k 2k 

migf* + ^(mng^) + S -^d kz (k z g^>) 



skzgo 



s<,> 



2k 
\m\ 2 6' 

2k n 



2k n 



9k z 9o 



(37) 



It is now easy to check that, to leading order in gradients, the Kubo-Martin-Schwinger 
(KMS) condition 



O k /T q< 



(3i 
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holds when n = uq. The first order solution (|3~2D with n = no does not obey the KMS 
condition however, except in the limit of a static wall v w — > 0. 

The solution for the scalar field Wigner function is much simpler. One can show |J 
that, to the order h in gradient expansion, there is no source for bosons, so that the 
bosonic Wigner function in the wall frame can be approximated by 



7T 

kty 



5(k - u^U + S(ko + uty)(l + /^) 
5(k - uty)(l + f^) + 5(k + u^ftf, 



(39) 



where /</,(fc z ) and f^,(—k z ) denote the occupation number of particles and antiparticles, 
respectively, such that in thermal equilibrium we get the standard Bose distribution func- 
tion 

1 1 

(40) 



?'1w{u<i > +v w k z )/T _ ^' 



(ffw(u<l>—Vwkz)/T _ l 



with 



k 2 + m\. Note that in thermal equilibrium the KMS condition, Aq 



e fco/T A^, is satisfied. 



5 Collisional sources 



5.1 Collisional source in the scalar equation 



Making use of equations (p2| ), (|28| ) and (p0|-p7|), one can show ||TT| that to first order in 
gradients the scalar field collision integral Cs is proportional to delta functions: 



(41) 



Because of the positive relative sign between and in the scalar Wigner function (|39D , 
the source arising from gives rise to a CP-even and spin independent contribution to 
the scalar kinetic equation (|22|) , which is thus of no relevance for baryogenesis. As a 
consequence, there is no CP-violating source for stops [|I^] in the collision term at first 
order in gradients. Our analysis is however based on a one-loop approximation for the self- 
energies. To make the study more complete, it would be desirable to perform a two-loop 
calculation of the source. 



5.2 Collisional source in the fermionic equation 

Before proceeding to calculate the source in the fermionic kinetic equation, we first prove 
that there has to be a nonvanishing source in the fermionic collision term. 



Since the equilibrium solution (P^| - PdT) for the Wigner function S ,<,> does not satisfy 
the KMS condition, the collision term does not vanish in the equilibrium defined by 
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Figure 4: Shown as a function of \m\/T is the expression |m| 2 Z^(|m|, m^), which appears in 
the collisional source (46) arising in the fermionic kinetic equation at one loop. For simplicity 
we have set T = 1. 



Ansatz (^). We shall now show that it is not possible to avoid this outcome by choosing 
n differently. 

First recall that for a static wall with v w = the equilibrium solution ( j32l - |36D does 



satisfy the KMS condition. Now try to modify the equilibrium solution for S l<,> as follows 

S< = 2iA(n + Sn), S> = -2iA(l - n - Sri), (42) 

where Sn is some scalar function of first order in gradients. The spectral function 

A = Aq + Ai (43) 

contains the spinor structure at leading and first order in derivatives, respectively, and 
the (^-function, 5{k 2 — \m\ 2 + s\m\ 2 9' /k ), as specified by equations fl3"7p . A\ contains 
derivatives which spoil the KMS condition when acting on rio(ko). We now seek 5n such 
that S > = —e /3h °S < be satisfied. Working to first order in gradients one arrives at the 
following condition 

A Sn = n Ail-Ain o . (44) 

where Ai acts as an operator on 1 and n , respectively. Since Aq and A\ have different 
spinor structure, there is no scalar function Sn that solves this equation. Hence there are 
indeed collisional sources that cannot be removed by a change of thermal equilibrium that 
is local in momentum. 



S 



In order to obtain the relevant collisional source in the kinetic equation for fermions, 
we need to multiply Eq. ( p7[ ) by — P s , take the real part of the spinorial trace and truncate 
at first order in gradients. 

To get the collisional source in the vector current continuity equation one should 
integrate the collision term over the momenta. It is not hard to see that the integrand 
is symmetric under the change of variables k z —>■ —k z and k' z — > —k' z , which then implies 
that 

r d 4 k „ 

W c *<°>=°- < 45) 

where C^(o) = —Ti[P s C^\ and f± denotes integration over the positive and negative fre- 
quencies, respectively. This means that there is no collisional source in the continuity 
equation for the vector current. 

The source in the first velocity moment equation, the Euler fluid equation, is however 
nonvanishing. One can show that, to leading order in v w , the source can be written in 
the following form 

[ d 4 k k z 2 s\m\ 2 6' 
2 Lj2^~ C ^ = wV 64^T /( ^' m '' m ^ (46) 

where the function Z/(|m|, m^) contains a complicated integral expression. To get a quan- 
titative estimate of the source, we perform numerical integration of and plot the result 
in figure |j. Note that for \m\ ^> T (T = 1 in the figure) the source is Boltzmann sup- 
pressed, while for |m| <C T it behaves as oc (\m\/T) 3 ^ 2 . The source peak shifts towards the 
infrared as m$/\m\ increases, and the magnitude drops drastically when one approaches 
the mass threshold = 2\m\, as requested for scalar particle decay and inverse decay 
processes. 



5.3 Collisional sources in the relaxation time approximation 

The methods used in literature for computation of sources from the collision term in scalar 



and fermionic equations |T2| , [L3|] can be in many cases rephrased as the relaxation time 
approximation: 

Cil>si{o) ~ ~ r» {fsi — fsio) > (47) 

where / S j denotes the true particle density of spin s and flavour i, f si0 is the equilibrium 
particle density (for a moving wall), and is the relevant relaxation rate, which is usually 
assumed to be given by the elastic scattering rate. In Paper II we have argued that, 
when the equilibrium distribution function is taken to be f s iQ = 1/ ^ UJai±+Vwkz ^ T + 1), 
the collision term (|47| ) contains no CP-violating source of relevance for baryogenesis. 
When one however takes (incorrectly) for the equilibrium distribution function / s jQ — > 
l/[ e (wo+vv)kz)/T _|_ -q^ implies the following form for the CP- violating contribution 
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to the collisional source 



1 i' L si — 



-r, : 



— n: 



r 



167T 2 



Jo(|M,|/T), 



(4£ 



where n° i± 



J ± (d A k/ (27r) 4 )(yfg and Mj and denote the masses and CP- violating phases 
of fermions. The detailed form of the function j7o is discussed in Paper II and it is not of 
importance for the purpose of this talk. 

Apart from parametric dependencies issues, the most worrisome feature about the 
relaxation time approximation interpreted such that it contains the source ([48), is that 
it yields a nonvanishing source (f4*8|) in the vector continuity equation, which we found to 
vanish based on equation (^). Moreover, the source does not vanish in the static limit 
v w — > 0, and hence it is clearly unphysical with no apparent relevance for baryogenesis. 
We finally note that the CP-violating sources 



discussed here, when applied 
to charginos in the MSSM, show the parametric dependence |Mj| 2 ©^ oc {h\h 2 )' on the 
Higgs field vev's hi and h 2 . This is to be compared with Refs. |TJL fTSI, where a source 



proportional to h\h! 2 — h 2 h[ was found and claimed to be important for baryogenesis. 



6 Conclusions 

We have studied collisional sources appearing at first order in gradients of a spatially 
varying mass for the model of chiral fermions interacting with a scalar field via a standard 
Yukawa interaction. This model is relevant for baryogenesis from fermions interacting 
with the Higgs condensate on growing bubbles at a strongly first order electroweak phase 
transition. The self-energies have been approximated by one-loop expressions. We have 
argued that, at first order in gradients, there is no collisional source in the scalar kinetic 
equation, indicating that baryogenesis sourced by scalar particles is highly suppressed. 
We have then proven that there is a CP-violating source in the fermionic equation and 
performed a quantitative analysis of the source. This source, together with the one from 
the semiclassical force, comprise the relevant sources for baryogenesis at the electroweak 
scale. Finally, we have argued that the collision term, when modeled in the relaxation time 
approximation, contains no CP-violating sources. In order to perform a full quantitative 
assessment of collisional sources, a two-loop analysis of self-energies is required, which is 
a work in progress. 
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